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Dstract 


‘Presented are experimental data obtained from the Ex- 
jlorer VIII satellite of five parameters pertinent to the prob- 
of the interaction of space vehicles with an ionized 
itmosphere. The five parameters are: photoemission current 
Jue to electrons emitted from the satellite surfaces as a result 
xf solar radiation, electron) and positive ion currents due to 
the diffusion of charged particles from the medium to the 
spacecraft, the potential of the vehicle relative to the medium, 
ind the ambient electron temperature. Included in the ex- 
perimental data is the aspect dependence of the photoemis- 
sion and diffusion currents. The observations then are used 
fo postulate certain characteristics of the satellite’s plasma 


sheath. 


Introduction 


_ This report is concerned with one type of interaction 
xetween a space vehicle and an ionized atmosphere. 
Specifically, it deals with the equilibrium potential of 
he Explorer VIII satellite and the current exchange 
xetween this spacecraft and the ionosphere. The ex- 
yerimental data herein presented are fundamental 
o the characteristics of the plasma sheath surrounding 
1 space vehicle. The report deals, then, with an inter- 
wtion localized to the proximity of the body. The 
lata should be important to the evaluation of other 
ypes of interactions that have been postulated since 
ll disturbances must begin right at the body itself. 
Both the characteristics of the ionosphere and of the 

pace vehicle contribute to their mutual interaction. 
‘or an orbiting satellite, the sheath properties are 
letermined by the following factors: 
\. Parameters of the undisturbed ionosphere 

(1) The ambient electron temperature (7’.); 

(2) the random electron current density (/<); 

(3) the random ion current density (J+); and 

(4) the magnetic field (B). 
3. Factors due to presence of the satellite 

(1) Radio-frequency fields used for telemetry 

transmissions; 
(2) conductivity of the surface ; 
1 Presented at the American Astronautical Society Meeting 
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(3) photoemission current density (J,) due to 
solar radiation; and 
(4) the satellite motion. 
For a conducting body at rest where rf and magnetic 
fields and solar radiation may be neglected, the equi- 
librium potential is given by: 


oo = ny (1) 


—kT, 
€ 
where k is Boltzmann’s constant and e the electronic 
charge. Since J. > J+, ¢ will be negative resulting in 
a positive ion sheath with a thickness related to T, 

and the electron concentration, N, . 

Consider now the factors introduced by the presence 
of the satellite. It is known that antennas radiating 
rf fields can, by rectifying process, cause a larger 
negative potential than would be expected from Eq. 
(1). However, experimental observations made from 
rockets (Ref. 1) during rf silence have been com- 
pared with data: taken in the same flight during rf 
transmissions. These data show that for the amount 
of power (100 mw) and the frequency (108 mc) used 
by the Explorer VIII telemetry system the rf field ef- 
fect due to telemetry transmissions can be neglected. 

The satellite motion affects the sheath in the follow- 
ing manner. The velocity is much greater than that of 
the positive ions. Consequently, the random ion cur- 
rent is essentially incident over that portion of the 
satellite surface projected in the direction of motion. 
The electron current should also be a maximum at 
the forward surface although the effect will be smaller 
because of the higher thermal velocities of the elec- 
trons. 

The electron current is also affected by the mag- 
netic field, as Beard and Johnson (Ref. 2) have de- 
scribed. The motion of the satellite with velocity V 
through the magnetic field produces an induced 
potential that is a function of position on the satellite 
surface, 


& = g + (V X B)-d, (2) 


where ¢) would be the satellite potential with no mag- 
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netic field, and d is the vector distance of any point 
on the surface from the satellite center. A satellite 
potential of do will be measured at all points which 
lie on a plane passing through the satellite center 
perpendicular to V X B. All other points will either 
be more positive or more negative than ¢o depending 
on which side of this plane they are situated. The 
electron current would be expected to be a maximum 
where ¢ is most positive which would be near the point 
corresponding to the direction V X B. 

The photoemission current density, Jp, tends to 
make ¢ more positive. At low electron concentrations, 
it is possible that photoemission can predominate 
resulting in a positive ¢9 and a sheath containing 
electrons. 

If the effects of the magnetic field and satellite 
velocity on the electron current can be neglected, 
then the following expression for the satellite potential 
is valid: 

kT 


= 2 Il 
do . n 


aS 


' (3) 
[148+ [7,48 


where the integration is over the satellite surface S. 
To take into account the variations in electron current, 
it is necessary to go back to the fundamental equation 
for current balance: 


[4.48 = [ 4,08 + | s,as. (4) 


This report will present experimental values ob- 
tained from the Explorer VIII satellite of the electron 
diffusion current, 7., the ion diffusion current 7, , 
and the photoemission current, 7, , as a function of the 
orientation of those points relative to the velocity, 
solar and magnetic field vectors. Also presented are 
measured values of ¢@ and 7’,. These experimental 
observations then are used to postulate a qualitative 
model of the plasma sheath and a quantitative mode] 
of the current exchange between the satellite and the 
medium. 


The Explorer VIII Satellite 


The Explorer VIII Satellite was launched from 
Cape Canaveral on 3 November 1960, onto an orbit 
with a 50° inclination, a perigee of 425 kilometers and 
an apogee of 2300 kilometers. It had a planned active 
life of two months. The primary mission of this Iono- 
sphere Direct Measurements Satellite, which is only 
incidental to this report, was the in-situ measurement of 
electron density and temperature and positive ion 
concentration and mass. 

A photograph of the satellite, highlighting the 
characteristics pertinent to this discussion is presented 
in Figure 1. The aluminum shell consists of two trun- 
cated cones joined at the equator by a cylinder. It is 
30 inches in diameter at the equator and 30 inches 
high. Non-conductive thermal coatings are located 
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Fic. 1. Ionosphere direct measurements satellite (E& 
plorer VIII). 


on both cones in a pattern conducive to the mail 
tenance of an equipotential surface. Retracted ai 
two ten-foot wires which served as a shortened dipa 
for an rf impedance experiment designed to measu 
N.. A combined solar-horizon seeker provided suj 
porting information on the satellite orientation. 

Data from only four of the satellite’s many senso: 

are considered in this presentation. They are: 

(1) an ion current monitor, responsive only to tk 
incoming positive ions; 

(2) an electron current monitor, responsive to tk 
sum of the incoming electron current ary 
outgoing photoemission current; 

(3) a total current monitor, responsive to all thre 
types of current; : 

and 

(4) an electron temperature probe, responsive on} 
to electron current and which measures 
ANC aes, 

The locations of these four sensors relative to th 
aspect sensor are shown in Figure 2. All but the electro 
temperature probe are centered on the equator. Th 
latter sensor is positioned near the forward end of tk 
spin axis as shown in Figure 1. 

In order to evaluate in detail the simultaneoy 
effect of all the sheath-influencing factors listed in th 
introduction, a single set of conditions was chose 
where favorable satellite orientation permits cled 
delineation of the dependency of the ion, electron, an 


| 
| 
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2A RETARDING POTENTIAL PROBE 
2B ELECTRON TEMPERATURE PROBE 
3A ELECTRON TEMPERATURE PROBE 
3B TOTAL CURRENT MONITOR 


. 2. Sensor locations (Explorer VIII Satellite) 


photoemission currents on the location of the velocity, 
solar and magnetic field vectors. Consequently, all 
of the reported data will be for the orientation described 
by Figure 3. Of major importance is that the velocity 
ind solar vectors are separated by 159°. This permits 
separate scanning of the solar and velocity dependent 
quantities as the satellite spins. The solar elevation 
mgle was 33° on the upper cone and the velocity 
vector depression angle was 15° on the lower cone. 
All the data were acquired within two minutes of 
22:57 UT on 27 November 1960, at which time the 
satellite altitude was 1000 km and its geographic 
soordinates 33°N and 84°W. The spin rate was 21.4 
rpm at this time. 


Measured Values of Positive lon Current 


The positive ion current flowing from the ionosphere 
(0 the satellite was monitored by a sensor shown 
schematically in Figure 4. The sensor is constructed 
n planar geometry and contains three parallel elec- 
rodes. The outermost grid is flush with and electrically 
sonnected to the satellite skin. The inner grid is biased 
iegatively to suppress photoemission from the collector 
ind to remove incoming electron current from the 
neasured collector current. This collector current 7 
s related to (7,)s, the current incident at the skin, 


ay 
14 = a+(t+)s, (5) 
vhere a, is the combined electrical transparency of 


he two grids for positive ions. 
Plotted in Figure 5 is the experimental collector 


MOTION 


Fig. 3. Satellite orientation (22:57 UT, 27 November 1960) 


SATELLITE SKIN 


TO TELEMETER 


Fic. 4. Ion current monitor (Explorer VIII Satellite) 


current as a function of the azimuth angle of the sensor 
relative to the velocity and solar vectors. The absence 
of a current when the sensor is pointed at the sun is 
proof that suppression of photoemission from the 
collector has been successfully accomplished. The 
behavior of 74 relative to the velocity vector is in good 
agreement with that predicted by Whipple (Ref. 3) 
in the following general equation which takes into 
account all values of the satellite-to-ion velocity ratio: 


i= aNeA lv cos 6 (3 + * orf r) + aoe i) (6) 


2 2/r 
where 
_ V cosé / ge 
a pe 
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and where A is the area of the collector, N the particle 
concentration, a the most probable thermal velocity 
of the particle, and 6 the angle between the sensor 
and the velocity vector. This expression has been 
derived for plane geometry and is valid for either 
ions or electrons if the appropriate values and signs 
for the symbols are observed. For angles less than 


SOLAR AZIMUTH ANGLE (DEGREES) 


about 45° in Figure 5 the observed currents are fitt 
by the reduced equation: 


i, = a4ANxeV cos 8. ( 


Using an inflight calibration of the combined electric 
transparency of the two grids (92% ), a value for ) 
of 1.3 X 10‘/em* is computed. This is consistent wi 
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electron concentrations obtained from the rf im- 
edance probe experiment also carried on Explorer 
II. Ionosonde data taken at this geographical 
osition and at this time yield an electron density of 
7X 10° el/cm’ at an altitude of about 300 km. The 
D0sitive ion concentration measured from the satellite 
erees with this value, assuming an ionospheric model 
ith a neutral gas scale height of 60 km, diffusive 
equilibrium, and a predominant 0~ constituent. 
One important qualitative conclusion related to 
heath characteristics can be made from Figure 5. 
‘he absence of a positive ion current on the side away 
from the velocity vector is definite experimental evi- 
dence for an electron sheath immediately adjoining 
the vehicle at this location. 
In Fig. 6 are plotted values of ion current density 
computed from 7,/A as a function of the total angle 
relative to the velocity vector. The solid line is a 
theoretical curve based on Equation 6, with ¢ = 0. 
The agreement between the observed currents and the 
BP istion is evident. However, in order to get agree- 
ment for angles greater than 45 degrees, a value of 
.6 X 10° cm/sec had to be used in the equation for 
The most probable ion thermal velocity. This value 
corresponds to an anomalously high ion temperature 
if only 0+ is assumed, but is reasonable if there is a 
significant amount of H+, implying that at 1000 km, 
the satellite is near the transition region between the 
upper ionosphere and the protonosphere. However, the 
most probable explanation for the fact that the meas- 
ured ion current for 6745° is higher than that predicted 
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Fig. 7. Electron current monitor (Explorer VIII Satellite) 


by Equation 6 for a reasonable temperature is that 
the negative satellite attracts ions that would otherwise 
not have been collected. 


Measured Values of Electron Current 


A sensor responsive to the sum of the incoming 
electron current and the outgoing photoemission 
current is illustrated schematically in Figure 7. It is 
mechanically identical to the ion current monitor 
but differs electrically in that the inner grid is biased 
positively rather than negatively. The positive bias 
serves to remove the incoming ion current from the 
measured collector current. The collector current is 
given by: 


1 = SF De == PACA =f Diners ; (8) 


where a, and a, are the respective grid transparencies. 
Plotted in Figure 8 is the experimental collector 
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Fic. 8. Electron current as a function of aspect (Explorer VIII Satellite) 
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current as a function of aspect. The current has its 
maximum positive value (photoemission) when the 
azimuth angle of the sensor relative to the sun is 
zero. As postulated in the introduction, the incoming 
electron current should be a function of the sensor 
orientation relative to both the velocity and the mag- 
netic field vectors. In this case, the maximum electron 
current on the shaded side is observed when the angle 
between the sensor and the cross-product of the velocity 
and magnetic fields vectors (V X B) approaches a 
minimum. This is as expected because at this time 
during a spin period the surface near this sensor loca- 
tion would be at its most positive point. 


Measured Values of Total Current 


The total current to the satellite was measured at 
its equator by the sensor illustrated in Figure 9. It 
consists simply of a collector flush with and insulated 
from the satellite skin. This collector current is given by 
= (as a Ce + Ge)e : (9) 
It represents then, the sum of the effects of Figures 5 
and 8, except that the overall amplitude of each 
component is larger since no grid transparencies are 
involved. 

The total current at the point of measurement is 
plotted as a function of aspect in Fig. 10. It is in 
definite agreement with the predominant features of 


i= 
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Fig. 9. Total current monitor (Explorer VIII Satellite) 
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Fie. 10. Total sheath current as a function of aspect (Explorer VIII Satellite) 
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-45 


135 


the graphs in Figures 5 and 8. Specifically, the curren 
peaks in the positive direction when the solar angle 1 
zero because of photoemission. The current peaks 
again in the positive direction for a velocity azimutl 
angle near zero. The peak here is due to the “ram 1 
effect”? upon which the positive ion current is largely 
dependent. The displacement from the zero angle point 
is due to the influence of the electron current in this 
region. Finally, a peak current in the direction of ary 
incoming electron current occurs when the angle 


velocity and magnetic field vectors (V X B) approaches 
a minimum. a | 
The satellite spin permits examination of the separate 


by comparison of Figure 10 with Figures 5 and 8 ta 
assign an electrical transparency to the grids used i 
the ion and electron current monitors. This in-flight 
calibration shows that the combined electrical trans- 
parency for positive ions of two parallel grids is ap- 
proximately that of the optical transparency. On the 
other hand, the electrical transparency for electrons is 
only about 30 percent of the optical transparency: 
This is an important observation in the evaluation off 
electron temperature data which are discussed in @ 
later section. 

A significant quantitative value derived from Fig. 10) 
is the measured value of the photoemission current 
density (5 X 10° amps/cm’) taken right at the 
vehicle surface for a minimum solar angle. This value 
can be compared with the random electron curren 
density computed from ionospheric models to predict 
an approximate altitude where the spacecraft potential 
can become positive. For most ionospheric models: 
this should occur at about 4000 km, not too far above 
the apogee altitude of Explorer VIII. 
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easured Values of Electron Temperature and 
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' The sensor illustrated schematically in Figure 11 
as used on a time-sharing basis to alternately monitor 
he incoming electron current at the sensor location 
and to obtain in the other half of its duty cycle electron 
temperature and equilibrium potential. It consists of 
two electrodes, a grid, and a collector. The collector is 
Bpeced positively to remove photoemission and_in- 
“coming ion current from the measured collector current. 
During one phase the grid is kept at the skin potential, 
thus permitting a measurement of the incoming electron 
current as a function of aspect. During the other time 
phase the grid potential relative to the satellite skin is 
‘varied from —1.2 to +8 volts to obtain ¢ and T,. 
The period of the grid voltage sweep (0.2 seconds) is 
‘kept small so that the collector current can ke studied 
as a function of the grid voltage for small changes of 
orientation. At the satellite spin rate, it was possible 
to obtain a volt-ampere curve for a change in orienta- 
tion of about 25°. 
__ Presented in Fig. 12 is a typical volt-ampere curve 
taken during the phase when the potential of the grid 
relative to the satellite skin was varied in accordance 
with the waveform illustrated in Figure 11. The 
electron current is plotted on a logarithmic scale. 
The shape of the curve is in good agreement with 
Langmuir probe theory. Two distinct slopes in the 
regions where the grid is below and above the plasma 
potential are apparent. The slope of the curve when 
‘the grid potential is negative is a measure of the 
electron temperature. The satellite potential is generally 
taken as the negative value of either the point where 
the curve departs from this slope or the point of in- 
tersection of the two slopes. For the illustrated curve, 
an electron temperature of 1800° + 300° K and a 
satellite potential between 0 and —0.15 volts is ob- 
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Fig. 11. Electron temperature probe (Explorer VIII Satellite) 
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Fic. 12. Typical volt-ampere curve for electron tempera- 
ture probe (Explorer VIII Satellite). 


tained. The error quoted in the temperature is due to 
the limited resolution. 

Not presented are additional volt-ampere curves 
taken at this time for different sensor orientations. 
These other curves show that as the satellite spins, 
the same values of electron temperature are obtained. 
The major difference is the current read at the plasma 
potential. These values are consistent with the electron 
current behavior when the grid is maintained at zero 
volts as shown in Fig. 13. 

The large scatter of points in Fig. 13 compared to 
all the previous graphs occurs because an electrometer 
sensitivity range is used where the input voltage to 
the telemetry system is small. The scatter, therefore, 
is due to random fluctuations of the telemetry sub- 
carrier at low signal levels. Despite this difficulty, it is 
apparent that the electron current peaks when the 
velocity vector azimuth angle is zero. From this ob- 
servation, together with an examination of the relative 
values of V and ¢ in Kq. 6, it can be concluded that the 
potential of the upper cone relative to the medium 
must be either close to zero as indicated above or can 
even become positive as the satellite spins. 

Figure 14 describes the orientation of the satellite 
with respect to the magnetic and velocity vectors and 
demonstrates that this implication is verified. The 
sensor on the upper cone always remains on the “‘posi- 
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Fie. 14. Orientation of satellite with respect to magnetic 
and velocity vectors (Explorer VIII Satellite). 


tive” half of the satellite, that is, the side that is 
positive with respect to ¢o. Moreover, the change in 
the distance from the ¢o plane is only 24 em correspond- 
ing to a change in potential of about 0.04 volts. In 
contrast, the potential on the equator changes by at 
least 0.14 volts and becomes more negative than ¢p . 
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Fic. 15. Qualitative satellite sheath model postulated from | 
experimental data. 


(Values of 7.4 X 10° em/sec for V and of 0.30 gauss} 
for B were used in the above computation.) Note also) 
that the electron current in Figure 13 is consistently) 
larger than that to the sensor on the equator, as one} 
would expect. . 
It remains to be pointed out that these conclusions| 
regarding the potential are necessary but not sufficient) 
to explain Figure 13. The observed current modulation | 
requires a lower electron temperature than the 1800° 
discussed in the preceding paragraphs. As suggested. 
previously, there is evidence that the electrical trans- 
parency of the grids to electrons is a function of the 
electron energy, the fast electrons getting through more. 
readily. This would tend to enhance any current modu- 
lation present, and it is suggested that this is occurring 
here. A changing grid transparency may also affect 
the determination of temperature. A thorough evalua- 
tion of the grid transparency effect on the electron 
temperature measurement is not yet completed. 


Sheath Characteristics Postulated from 
Experimental Observations 


_ It is concluded that the experimental data from the 
various sensors are mutually consistent. These data 
can be used to postulate a qualitative model of the 
sheath characteristics and a quantitative model of 
the current exchange between the satellite and the 
medium. 
_ First, the value of the satellite potential on the upper 
cone is between 0 and —0.15 volts. The satellite surface 
‘is not equi-potential. The motion of the satellite 
‘through the magnetic field causes an induced potential 
so that the potential is very close to zero or positive 
“in the V X B direction and correspondingly more 
negative at the other end. 
_ The potential distribution on the surface has a 
marked effect on the observed electron current, with 
most electrons incident at the more positive end. 
It should be noted that there was a net electron current 
_ to the skin even when the current flow was perpendicu- 
lar to the magnetic field. A qualitative model of the 
sheath surrounding the satellite is illustrated in Fig. 
B15, showing also the effects of solar and velocity vector 
orientations on the various current exchanges. The 
_ electron current shows a ram effect due to the satellite 
motion if the frontal surface is favorably situated with 
respect to the more positive end of the satellite. 
_ The positive ion current density is a function only 
_ of the angle between the surface normal and the satellite 
velocity vector. The absence of this current behind 
_ the satellite indicates that there is an electron sheath 
or wake adjoining the rear surface. If it is assumed 
that the wake is in the form of a cone, its size may be 
estimated from the ratio of satellite to ion velocities 
and the satellite diameter. In this case the cone has a 
half-angle of about 25° and extends back a distance of 
about one satellite radius. A positive ion sheath sur- 
rounds the front of the satellite and should envelope 
both the satellite and the electron wake. The thickness 
of the ion sheath would be comparable to one Debye 
length which is computed to be 2.5 cm. 

Finally, there is a photoemission current on the 
satellite surfaces facing the sun. This current is im- 
portant in considering the current balance to the 
satellite and thus the equilibrium potential. However, 


photoemission does not appear to disturb the sheath 
adjacent to the emitting surface. 

The sheath model illustrated in Figure 15 is for an 
altitude of 1000 km. The measured positive ion density 
is 1.3 X 10°/em*. A neutral gas scale height of 60 
km corresponding to a temperature of about 1100° K 
is indicated by comparison of this N* with ionosonde 
data. An electron temperature of 1800° + 300° K 
was measured by the electron temperature probe. 
However, this value could be influenced by the effect 
on the probe’s grid transparency of a changing electron 
retarding potential. The percentage modulation of the 
electron: current by the satellite velocity indicates a 
somewhat lower value for 7’,. The behavior of the 
positive ion current at the satellite sides suggest the 
possibility that the transition region between the 
ionosphere and the protonosphere is near 1000 km. 
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Furthering Basic Biological Knowledge 
Through Space Research 


June Lee Biedler? 


Abstract 


The advent of the ‘“‘space age”’ has resulted in both the op- 
portunity and the necessity for extended research in various 
areas of science. That biological science is gathering impetus 
in certain directions, some of them new, is illustrated by the 
results of a survey taken for the purpose of this report. There 
is a manifest need for an interdisciplinary approach to space 
research whereby those working in the biological sciences will 
need to further their acquaintance with engineering possi- 
bilities and engineers with the requirements of biologists. 


Introduction 


It seems timely to assess the possible contributions of 
space exploration to the furtherance of basic biological 
knowledge—knowledge which is derived from the ana- 
lytical or observational study of life processes, be they 
in unicellular organisms such as bacteria or certain of 
the algae, or in higher plant or animal forms such as 
tomato plants or mice. While not excluding man as a 
worthy representative of the biological world, in general 
the medical or physiological, social or psychological 
considerations of man-in-space are not pertinent to this 
discussion. Many of the present studies and projected 
plans are directly related to that goal—the survival and 
sustenance of man in the space vehicle or in the extra- 
terrestrial environment. But it is the other side of the 
coin which concerns us here—the fundamental knowl- 
edge which is now increasing and which can be extended 
by the pursuit of this objective. 


Space as a Biological Laboratory 


Just what are some of the physical factors unique to 
the space environment that may make possible an ex- 
tended understanding of certain biological mechanisms, 
or may even lead to entirely new concepts? There are 
kinds and intensities of cosmic and solar radiation, 
various planetary environments, high vacuum to the 
order of 10~'8, zero-gravity, and so forth. The possibili- 
ties of space research have excited those few biologists 
who have followed the development of space vehicles, 
not because such research will necessarily bring a simple 
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solution or sudden enlightenment to persistent bio- 
logical questions, but because there is provided a means 
of observing terrestrial life mechanisms away from a | 
tain constant and inescapable influences in the earth | 
environment, as well as a possibility of discovering j 
mechanisms which have originated entirely outside of 
these influences. 


Implications of Zero-gravity 

The physical factors involved are well known to. 
engineers and other physical scientists. Possibly the 
most important of these, at least at this incipient phase | 
of the ‘“‘space age,’ is zero-gravity, the reduction or | 
complete absence of the gravitational force. We are able | 
to reproduce the converse of this state, an increase in 
gravitational pull, by means of the centrifuge. However, | 
we cannot with validity extrapolate back to the zero- | 
gravity state. Nor can we with complete surety and | 
accuracy simulate this state, although attempts at an | 
approximation to zero-gravity are being made and are | 
worthwhile. Of the greatest importance at present is | 
the effect on man of this entirely strange factor. But how | 
can we begin to understand the possible or even the ob- | 
served effects upon the highly organized, complex | 
human being before we are able to know if or how | 
certain readily observable biological phenomena are 
altered. . 

What effect has zero-gravity on cellular division? For — 
example, it is known that in a given nutritional milieu : 
at a controlled pH, at a constant 37° C temperature, a 
certain strain of human cancer cells grown in glass cul- 
ture flasks will divide with an average generation time 
of 22 hours. Will a change from 1g to Og alter this di- 
vision time? Some biologists would say yes, probably; 
while others would say no, there could not be the 
slightest effect. The point is, this question cannot be 
answered on earth. There is a strain of bacteria, 
Escherichia coli, commonly present in our intestinal 
flora, which will, under optimal conditions, divide ap- 
proximately every 30 minutes. Would this division 
time be affected? Certainly a simple question and a 
relatively simple one to answer from the point of view 
of experimental design for a space vehicle; but has this 
question been as yet satisfactorily answered? 

What is the effect of zero-G upon embryological de- 


velopment? There are some classical experiments which 
tt would be well to reassess under new conditions. For 
xample, the embryologist may turn the fertilized frog 
g¢ upside down to see if the change in polarity of the 
early embryo influences the development of the tadpole. 
t turns out that during a certain critical period at least 
it does—the cleavage of the egg and the tadpole, if 
ormed at all, are abnormal. Exactly why is not known. 
hat would be the effect of the elimination of the con- 
stant gravitational pull? It is a fascinating question, 
and the answer may well illuminate this little under- 
stood problem. 


Present State of Biological Science 


_ In the foregoing we have discussed the biologists’ 
interest in extending knowledge of biological processes. 
It might be instructive to mention the state of the art 
of biological science as compared with that of the 
physical sciences, since this bears on later remarks con- 
cerned not so much with the relationship between the 
biological and physical sciences as with the communi- 
ation between biologists and the physical and engineer- 
‘ing scientists. 

- Perhaps it is not widely realized how great is the 
abundance of uncertainties and unknowns in the un- 
derstanding of life processes. In biology there are few 
general laws, whereas physical science has advanced to 
the state where there is a great body of laws which may 
be expressed in mathematical terms and by which phe- 
nomena may be predicted with complete accuracy. A 
well-known characteristic of biological phenomena is 
their variability. Lack of perception of biological un- 
knowns on the part of engineers and those planning 
space missions may result in insufficient lead time for 
the gathering of biological data upon which orderly pro- 
gression of man-in-space missions depends. Credit is 
given to biologists and medical scientists for possessing 
information or understanding, credit which unfortu- 
nately is as yet so often undeserved. This is the reason 
why there can be few ready-made answers to problems 
posed by even our next steps in space flight. Compound- 
ing the problem may also be political pressures, which 
mitigate against the freedom always necessary to the 
pursuit of basic research. 

There needs to be discernment. and definition of the 
objectives of research in space, and perhaps the realiza- 
tion that the accumulation and dissemination of bio- 
logic and other scientific information is of tremendous 
import to the prestige of this nation, as it is to any 
nation. Despite the potential high interest of the biolo- 
gist in the utilization of space exploration to accumulate 
knowledge on basic life processes including the origin of 
life itself, there is a mere handful of biologists in the 
United States who are thinking in terms of this aspect 
of our future. Responsible in part for this is the lack of 
realization on the part of biologists, immersed in their 
own research problems, of these new parameters, the 
opportunities afforded by the non-earth environment 


for the extension of understanding. In addition, from an 
ethical point of view, the biologist may be deterred from 
active participation by the enormous sums of money 
involved in space vehicle development. And this is an 
attitude undoubtedly shared by biologists and non- 
biologists alike. Here is an area in which the engineer 
can be helpful by elucidating the capabilities of pro- 
jected space vehicles in order that biological experi- 
ments may be designed so as not to impose a significant 
additional expenditure. 


Necessity for Fundamental Research 


Despite certain present-day deficiencies in planning 
and outlook, as noted here, we may be encouraged and 
stimulated by the studies that are now being conducted. 
For this report a survey was made of biological pro- 
grams related to the space effort, and many of these are 
briefly described below. A number of investigators, from 
whom information was collected concerning experi- 
mental projects underway, replied quite candidly that 
their work preceded interest in space exploration or was 
in effect work of fundamental significance with implica- 
tions for space exploration only secondary. A large 
measure of credit must be given to NASA and to the 
Department of Defense who are supporting these 
studies. Most of these experimental programs are quite 
new and there are few results as yet. Of course there is 
another, rather more impelling reason why it is neces- 
sary to speak more of plans and purposes than of results: 
the fact that very little biological experimentation in 
space has as yet been possible. 

Research Projects 
Zero-gravity 

Dr. Charles Lyon at Dartmouth College feels that in 
the plant world the fast-growing herbaceous annuals 
might well be found to be the most efficient producers 
of oxygen as well as supplementary food sources. Since 
each crop must be started anew from seed, the effect of 
gravity-free space as well as other space capsule condi- 
tions upon seedlings must be checked. He approaches 
this problem from the point of view that since the 
growth of plant tissue is quite closely controlled by the 
auxins, which are growth hormones, any gravitational 
effect upon auxin transport might play an important 
role. Lyon has demonstrated a marked effect of gravity 
on plant growth by the use of the clinostat, long used 
by botanists as an experimental means of obviating the 
effects of normal gravitational action. 

From Emory University, Professor Geoffrey Bourne 
has submitted a project for the study of the weightless 
state on muscles in tank immersion experiments as well 
as actual zero-gravity on muscle fibres and other cells 
as well. 

Dr. Cornelias Tobias and others at the University of 
California are concerned with the effects of the gravity- 
free state. For example, they are interested in a beetle, 
Tribolium confusum, the adult and larva of which infect 
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flour and grains, since it may, in a space flight situation, 
give some information as to the effect of weightlessness 
upon differentiation. In general, they wish to evaluate 
those aspects of physical chemistry which will influence 
the metabolism and proliferation of cells in a state of 
partial or complete weightlessness. Some of the parame- 
ters of living cells that may be altered are aerobic and 
anaerobic metabolism, induction of mutations, rate of 
proliferation, differentiation, and morphology. 

Dr. Stephen W. Gray and his group at Emory Uni- 
versity replied that the nature of their overall interest 
in space research was to gain basic understanding of the 
role of gravity in determining growth and form of tis- 
sues and cells. In conjunction with Dr. Betty F. 
Edwards he has published results of the effect of 
supranormal gravitational forces on cell division in 
wheat seedlings and on the swimming habits of tad- 
poles, and is studying the effects of these forces on chick 
heart muscle grown outside the body under laboratory 
conditions. He emphasized that the most important 
gravitational problem for immediate pursuit is the be- 
havior of cells and tissues in gravitational fields less 
than that of the earth. Whereas high gravity may be at- 
tained experimentally, sub-normal and zero-gravity 
cannot, and to extrapolate from known effects in high 
fields to expected effects in low fields is dangerous. Long 
exposures to zero-gravity may be necessary; it may be 
that growing tissues will be affected while mature tissues 
will not. 

Dr. Charles Wunder at State University of Iowa is 
primarily concerned with gravity as an environmental 
agent, which, although almost constant on earth, would 
have different intensities at other locations within our 
solar system. For several years he has been observing 
the growth of fruit fly larvae, mice, and hamsters during 
continual centrifugation, in order to elucidate the role 
of gravity in the control of growth. It has been found 
that continual exposure to high gravity, as simulated 
by centrifugation, can reduce the growth rate of fruit 
fly larvae. Preliminary results show that between the 
first and third months after birth, mice, in a continually 
operating centrifuge, exhibited growth inhibition at 
2 g’s and 5 g’s and did not survive at 7 and 10 g’s. 

De James C. linn, Jr., at the Aerospace Labora- 
tories, Missile Division, of North American Aviation is 
in the process of studying the role of constant accelera- 
tion in growth and development of higher plants, with 
the purpose of determining the effects of simulated zero 
to 2 g’s on growth and development. 

In 1958 Dr. H. J. Muller, an eminent. geneticist, 
offered a system of providing a pseudo-weightless state 
for human beings, patterned after the clinostat, which 
has been mentioned as being used for obviating gravita- 
tional effects on plants. Several others simultaneously 
came forth with this basic idea, the nullification of 
gravity-orientation cues. Dr. Raphael B. Levine at 
Lockheed has developed a null-gravity simulator for 
human beings. But to reiterate, the gravity-free state 
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can only be simulated. Such methods will give somilj 
information of very practical value for man under 
earth-free conditions. However, fundamental cuestional 


concerning the influence of zero-gravity on biologic 
processes simply cannot be answered in this way. 


Radiation 


The interaction of living material and various types | 
of electro-magnetic and corpuscular radiation produces 
a broad spectrum of physical, chemical and biological 
effects which have been under study for a little more. 


than half a century. The discoveries of Becquerel, 
Réntgen, and the Curies opened up whole new areas in 


medicine and the natural sciences. Today the destruc- | 
tion of both malignant and normal tissue with the use | 
of X-rays and radioactive substances is a well-known | 
if not fully understood phenomenon. In present day 


cellular biology the uses of various forms and intensi- 
ties of radiation are manifold. 

One of the most fascinating of these is the production 
of mutations by the alteration of a finite part of the 
self-replicating, hereditary material of the cell, the | 
DNA. The observation of the specific effects, be they 
gene mutations, gross chromosomal alterations, or more 
extensive cellular damage, of a specific kind of radiation 


| 
| 


can provide valuable information as to the structure or 
function or composition of mammalian cells, bacteria, — 


plant cells, viruses, and so forth. Such studies are of 


basic importance, as well as of immediate practical 


value to man. 

Dr. Tobias and his group have for several years been 
studying the biological effects of heavily ionizing radia- 
tions as well as their pertinence to space flight. In addi- 
tion to his gravity studies, Dr. Bourne has initiated a 
project for the study of the effects of the types of radia- 
tion met by high altitude rockets and satellites on the 
muscular function of live frogs. Dr. Franklin Hutchin- 
son at Yale has been interested in the biological action 
of very heavy ions. The chief purpose has been to gain 
further understanding of the living cell and cellular con- 
stituents. However, since fast heavy ions are found in 
primary cosmic radiation, the importance of this work 
to space travel has become apparent. One aspect of 
these studies is the effect of these ions on chromosome 
structure detected by the production of chromosomal 
aberrations. 

Dr. A. Gib DeBusk of Florida State University has 
an extensive program for research in space genetics, 
which involves the exposure of suitable biological sys- 
tems to the Van Allen radiation belts and to solar 
cosmic radiation over a prolonged period of time. Very 
elegant genetic bioassay systems have been devised as 
an essential preliminary to man’s entrance into these 
regions. Dr. DeBusk is interested in answering the ques- 
tion as to whether primary cosmic radiation ean produce 
gene mutations. It is known that they cause cell death, 
but do they cause mutations, either in the molecular 
Sense or as gross chromosomal changes? Of timely 


interest is his report that on September 19, 1960 his bio- 
logical test system, using the bread mold, N eurospora, 
was exposed to the inner Van Allen radiation belt, and 
that this was the first so-called biological flight. The 
ample was recovered, having reached an altitude of 
ipproximately 1130 miles for 26 minutes, and is under- 
oing genetic and physiological tests in his laboratory. 


Algae and Other Plants in a Closed Ecological System 


_ The algae have come under close scrutiny as a possi- 
ble biological system for maintaining man within the 
closed ecological environment of the space vehicle. 
Although the purpose here is an eminently practical 
one, such studies lead to the accumulation of further 
; basic information concerning this large and diverse 
group of organisms, some of which are unicellular, some 
t of which are in aggregate form, as for example seaweeds. 
Moreover, great impetus has been given to the study of 
: BP otoeynithesis: of fundamental importance to our 
existence. 
_ The algae have shown the greatest promise for on one 
hand acting as a photosynthetic gas exchanger and on 
the other as a source of food for human beings. The 
presence of chlorophyll permits, in simple terms, 
thermodynamic conversion of light energy into chemi- 
3 ‘eal energy, whereby CO, and water are reduced to car- 
bohydrate with the liberation of oxygen. The initial 
product is sugar which is converted and stored as 
starch, upon polymerization. Other materials, such as 
proteins and vitamins, are formed indirectly. 

Dr. Robert Krauss at the University of Maryland 

“initiated two research programs concerned with certain 
aspects of the basic physiology of the algae. His group 
was the first to grow algae in a closed recycling system 
for 3 months with only inorganic nutrients. It was a 
critical test of the absorption rates of inorganic ele- 
ments, and it contributed to the present concept of the 
use of algae in the closed ecological system. He is inter- 
ested in the effect of the space environment at a cellular 
level, but points out that the studies are concerned pri- 
marily with biological phenomena of fundamental 
interest. 

Dr. Hiroshi Tamiya of the Tokugawa Institute of 
Tokyo has conducted a number of studies on the use of 
algae, primarily Chlorella and Scendesmus which belong 
to the group of green algae, as food or animal feed. In 
1942 several German scientists first suggested the large 
scale industrial production of algae as a food source. As 
Dr. Tamiya points out, the use of algae as food has been 
practiced in the Orient since ancient times, especially in 
Japan. Dr. Tamiya’s studies have not been concerned 
with the use of algae in the space vehicle but rather with 
basic growth and nutritional problems. Properly proc- 
essed Chlorella is non-toxic to human beings and has 
been found to be rich in proteins and vitamins as well 
as to contain carbohydrates, lipids, and minerals. Com- 
pared to animal protein its greatest drawback is its low 
content of one of the essential amino acids, methionine. 


Another interesting project is the growing of Chlorella 
with a species of Daphnia or water-fleas, by which it is 
hoped to devise an efficient system for the conversion of 
plant protein to animal protein. 

An extensive program of biological and engineering 
studies is underway at the Electric Boat Division of 
General Dynamics for the development of a photosyn- 
thetic gas exchange system for space vehicles. The 
algae are being utilized, at least at present, although the 
plant kingdom is being surveyed in detail. A particular 
strain of Chlorella, discovered in 1953, has been selected. 
It is a thermophilic strain, growing and photosynthesiz- 
ing optimally at 39 to 40° C. and with a tremendously 
increased growth rate. Heavy cell concentrations of 
this alga may be employed, since the requisite light 
sources for a maximal rate of photosynthesis can be of 
a very high intensity without inhibition due to light 
saturation. An additional project is the wide-scale 
screening of algae in order to find an even more efficient 
and practicable organism. 

Dr. Robert Tischer at Mississippi State College is 
doing basic research on the physics, chemistry, and mi- 
crobiology of closed ecological systems, in order to 
elucidate possible biochemical pathways from human 
waste products back to human food. Various kinds of 
microorganisms are under evaluation, for example cer- 
tain bacteria, molds, yeasts, and protozoa, as possible 
recycling agents for the provision of nutrient for algal 
cultures. 

At North American Aviation research is in progress 
to examine the possible role of the angiosperms, higher 
plants, in space operations. The stated purpose is to 
determine the optimal growth and development of cer- 
tain plants in reference to gas exchange rates and food 
production under conditions under which all environ- 
mental factors, except magnetism and high energy 
radiation, are controlled. 

Dr. Frank Salisbury and Dr. Ralph Baker at Colo- 
rado State University are interested in various aspects 
of plant growth under artificial environments. They feel 
that higher plants can be made to photosynthesize with 
greater efficiency, that they are generally more edible 
than the algae, that they might be useful in a water 
purification system, and that in any case their utiliza- 
tion is almost certain for future space exploration or 
planetary habitation. Two projects have been initiated. 
One is the response of certain plants and their photo- 
synthetic rate to high light intensities as might be found 
on the moon, in addition to other variables such as ele- 
vated temperature and a high CO, atmosphere. The 
second project is the study of plants grown under germ- 
free conditions, similar in concept to the animal studies 
which have been underway for many years. As is well 
known, associated with the intestinal tract of animals 
including human beings are various kinds of non- 
pathogenic bacteria. The roots of plants likewise exist 
in close association with a variety of bacteria, fungi, 
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and so forth. Their influence is considered significant, 
but has not been well elucidated. 


Biological Rhythms 


The subjection of men to alien environments wherein 
the familiar alternation of daylight and darkness are 
absent, where he is outside the earth’s magnetic field, 
and where he might not be able to have his approxi- 
mately 8 hours sleep every 24 hours, has led to a 
heightened interest in the many rhythmic processes in 
plants and animals which are, or seem to be, related to 
the 24 hour solar day. There has been a considerable 
amount of research in this area since the early *30s, 
chiefly in Germany and Holland as well as in this 
country and the United Kingdom. In broad perspective 
the problem is this: Are the many observable daily 
biological rhythms built into the organism as a sort of 
biological clock, or are they controlled by an external 
geophysical periodicity? The basic studies of such in- 
vestigators as Dr. Colin 8. Pittendrigh at Princeton and 
Dr. Frank A. Brown, Jr. at Northwestern University 
are of tremendous significance to space travel. In this 
area, likewise, Dr. Finn at North American is conduct- 
ing a study concerned with this question of endogeny 
vs. exogeny of biological rhythms in plants, in order to 
determine whether diurnal fluctuations in the earth’s 
magnetic field are responsible for the periodicity (but 
not the phasing) of daily biological rhythms. 


Other Experimental Areas 


There are many other areas of fundamental research 
which do not fall conveniently into the previously men- 
tioned categories. For example, Dr. Emmett W. 
Chappelle at RIAS is interested in the fundamental 
mechanisms of carbon monoxide fixation—that is to 
say, the incorporation of CO gas into a compound 
utilizable to a living organism. There are several micro- 
organisms which are known to be able to do this, among 
them the green alga Scendesmus. Carbon monoxide fixa- 
tion may be of practical importance in closed ecological 
systems. 

Dr. L. Joe Berry at Bryn Mawr College is concerned 
with the needed additional knowledge as to the response 
of human beings to high altitudes or lowered oxygen 
tension over a long period of time. It is postulated that 
animals carrying pathogenic organisms can become 
overtly diseased at lowered atmospheric pressures. The 
effect in mice of various combinations of altitude and 
oxygen concentration on tissue citric acid, carbohydrate 
lability, susceptibility to infection, and other physio- 
logical processes are being investigated. 

Extra-terrestrial Life and the Origin of Life 

This last category is of great significance to natural 
science as well as highly stimulating to man’s imagina- 
tion. But while there are a number of projected studies 


there is obviously little positive information with which 
to work. The difficulties are compounded by the fact 
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that we do not truly know what to look for. On earth 
all living material has a backbone of carbon, which is 


derived from CQ, and which is converted to organi , 


during the photosynthetic process. But did the atmos+ 
phere always contain CO,? And what happened before: 
the highly complex chlorophyll molecule was elabo~ 
rated? Then of course there is the all important ques+ 
tion of whether there could be a self-replicating systems 
based on an element other than carbon, and whether w 
would recognize it as a ‘‘living” system. 

There are at present several closely related actives 
programs. Since it is highly possible that under suitable: 
conditions, life might have arisen independently on 
other planets, recommendation has been made by thet 
International Committee on Contamination by Extra-- 
terrestrial Exploration that space vehicles be sterilized! 
in order to avoid contamination of the moon and other’ 
planets. In this way confusion of results from life detee- 
tion systems as well as possible harmful effects of the: 
organisms themselves in an alien environment can be} 
avoided. We know that the Russian moon probe was: 
reportedly sterilized. In the United States there are) 
rigorous sterilization procedures and active programs 
for their further development. This is in itself a difficult. 
task and necessitates the close collaboration of biologist. 
and engineer. Dr. Joshua Lederberg, the Nobel Prize 
winning microbiologist at Stanford University, has con- 
cerned himself with this problem, especially in relation 
to his plans for a cytochemical study of extra-terrestrial 
organisms by the use of a vidicon microscope for tele- 
metered observation. Elsewhere other detection sys- 
tems are being developed which involve, for example, 
the examination of planetary materials, or the analysis 
of planetary reflection spectra. 


Conclusion 


It was hoped that some suggested requirements, in 
engineering terms, for conducting biological experiments 
in space, could be reported. Though this question was 
posed specifically to investigators in the course of 
collecting information on present biological experi- 
ments, actually only two replies were volunteered. 
Perhaps the reason this was so has already been brought 
out and can well be reiterated here in summarizing. 

The development of a particular crossbreed of life 
scientist and engineer specifically for the design or per- 
formance of extraterrestrial biological experiments, as 
well as for manned space flight, may be undesirable in 
an age of needed specialization. However, until this 
interdisciplinary capability becomes the general case, 
rather than exhibiting itself in only a few especially 
talented and motivated individuals, we are under an 
imperative demand. The communication barriers be- 
tween the biologist and engineer, or generally between 
the life scientists on one side and on the other those 
designing space vehicles and planning space missions, 
must be broken down. 


The efforts of independent organizations such as the 
erican Astronautical Society; the support being 
iven by the Department of Defense to attracting the 
erest of biologists; the creation of the NASA Life 
ences Division and the drafting of its 10 year pro- 
ected program are all so guided. However, recognition 
the problem and educational effort exerted not only 
ough the public news media but also through the 
particular journals, or reviews, read by biologists and 
engineers may accelerate the kind of thinking which 
peeds up lead times in experimental concept and de- 
sign. The widespread dissemination of even preliminary 
results will stimulate both thought and action. Further- 
pore, international cooperation, the blending of our 
engineering capability with the very real and creative 
abilities of scientists of other countries, would be un- 
Measurably productive to biological and medical sci- 
ence, as well as contributory to understanding between 
peoples. 

Two specific points have been advanced in this report 
by which the engineers and the physical scientists could 
be guided in these space age endeavors. One is that bio- 
logical science is still in a comparatively early develop- 
mental stage. The other is that a clear grasp is needed 
of those areas of space flight which require biologists to 
be an integral part of the earliest planning team, so that 
the necessary lead time is provided for observation, ex- 
perimentation, and study, that space exploration can 
progress in an orderly, scientific fashion without grave 
and unnecessary stumbling blocks. 

In conclusion, it is suggested that man’s role as ob- 
server and monitor receive frequent and penetrating 
re-evaluation. From the descriptions given here of what 
are the very early space researches conceived by scien- 


tists, it is a short step to imagining the vast possibilities 
of future scientific experimentation in space or plane- 
tary environments. And it is hoped, from this brief 
exposition, that there is increased awareness of the very 
real advances that are being made right now as well as 
the tremendous potential that is offered for the further- 
ance of basic biological knowledge. 
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The Orbital Motion of Pellet Clouds 


Stanley Ross? 


Abstract 


The motion of a cloud of pellets is analyzed following their 
sudden release from a spinning container. The incremental 
velocities imparted by the spin motion are interpreted as time 
mtegrals of impulsive perturbation accelerations which modify 
the parameters of the container’s orbit. By transforming these 
yrbit element changes into linear displacements, the subse- 
yuent cloud shape may be studied analytically. Some exam- 
yles are included which illustrate the results for certain 
special cases. 
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Palo Alto, Calif. 


Symbols 

At the instant of pellet release: 

Symbol Definition 

o particle angular position in the container 

b spin velocity of the container 

p particle vectorial position, measured from some fixed 
point on the spin axis 

V= Xp 

a yaw angle of container, measured positively to starboard 

8B pitch angle of container, measured positively upward 


For the container’s orbit: 
Q,7,4,e,n,€, are orbital parameters, as defined in Ref. (1). 
v instantaneous true anomaly 


uU= v+te 
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h.= esnw 

l= e os w 

R, . W, “are unit perturbing accelerations, directed radially 
outward, circumferentially forward and orthogonally 
upwards, as defined in Ref. (1), page 402. 


ip time, measured from instant of explosion 
dz, dy, dz are rotations measured about orthogonal reference 
axes. 
fae 
na 
Introduction 


One means of creating an orbiting aggregate of 
individual particles is to suddenly burst a spinning 
shell packed with small pellets. The freed particles, 
each experiencing an additional velocity component 
imparted from the spin motion, separate to form a 
cloud which continues to circle the Earth. 

Suppose the angular velocity of the casing is specified 
by the vector ; then a particle which occupies a 
position p measured from some fixed point on the 
spin axis acquires an incremental velocity V = d X p. 
The motion of each pellet prior to the shell explosion is 
described by the orbital parameters common to the 
aggregate, e.g., Q, 7, a, h, 1, €. Removal of the container 
produces changes in the orbit of each particle which 
may be characterized by the familiar equations of 
variation of the elements. If the rupture is assumed to 
occur instantaneously, the ‘‘perturbing accelerations” 
may be treated as impulse functions whose time in- 
tegrals are merely the appropriate incremental velocity 
components. The orbital elements then change in 
accordance with the following relations:’ 


i= | oe | (v.W), 
na — e’ sin 2_)0 
(1) 
: r CO A 
in = — | (V-W), 
Ex - 5| 
2e sin v 2 
aa =| Seine | (Vv. Ry +[7*Y Zon Ao] (V-S), 
nvV/1 — elo 
hes | Sabet e COs | (V-R) pv — é sin 4 
na 0 na 0 


— 2 e@ A 
= = E — cos u sin | (V-S) 
i} 


na(1 — e?) 


re sin U COS w A 
— le a oe] (V-W), 
nav/1 — & tan ajo 


Hise: ee = e sin ‘ (V-R) re oe — é cos "| 
0 


na 0 na 


* Anticipating that shell orbit eccentricities will be small 
for many problems of practical interest, we have chosen the 
variables h, 1, « to specify the orbit in view of their analytic 
behavior near e = 0. That is, h = e sin w, 1 = e cos w, 
« = 7 — o (Ref. 1, pp. 406 and 421). 

“We assume that K, the ratio of incremental speed to 
orbital speed, is small. 
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POINT OF SHELL /@ 
EXPLOSION; u=u= 01 


Fia. 1. Resolution of element perturbations 


—., . 
+ Be sin wu sin of (V-S) | 


na?(1 


+/ re sin U COS w | (V-W), 
nar/1 — e? tan ito 


and | 


‘i = E ae avil—é (1 — V1 — &) cos of} 
na é 0 


-(V-R) Dl eee sin v (1 +) d-v1— | 


h r sin wu tan 5 s 
CV ata eee ) 

Insite, Go Ros: W and are unit vectors in th 
radial, circumferential, and orthogonal direction 
(Ref. 1, p. 402), while the other quantities on the righ 
sides of Eq. (1) (the symbols are those of Ref. 1 
are evaluated at the moment of explosion. 

Since we are examining two-body motion about | 
spherical Earth, there is no physically preferred diree 
tion in space. Therefore, any coordinate system ma: 
be chosen which will simplify the analysis. If a referene 
equator is established orthogonal to the shell’s orbi 
and passing through the position of the explosion (se 
Nig. 1), then dQ, dz, and dw may be transformed int 
small rotations dx, dy, and dz about the axes shown: 


dx = —cos udQ + sin u di, 
dy = sin udQ + cos u di, 
(2 
and 
dz = dw. 


Also, since 
_ _a(l —&) 
1 + e cos(u — w)’ 


1e radial displacement at any time t may be expressed 


_ (* an ar 
or = =) da + (2) de + eh dao. (3) 


Finally, by adopting the formula in Ref. 1, p. 171, 
nd from the definitions of h and | (n.b., uw = 0), it 
can be shown that 


n(t—t) =u—2[lsnu—h(cosu—1)]}+--- (4) 


where % is the time of explosion. For T = t — t, 
therefore, the angular gain in particle position along 
the orbit is given by 


dz+ du = dw 


2{dl sin u — dh(cos u — 1)| — et T da 
as a) 
1 — 2(1 cos wu + h sin u) 


= 
Thus, at any time after the release of the pellets, the 


positional change of each particle may be closely 
approximated by the following expressions: 


: CA, Ee 
dS = ridz + dul, 


‘> 


(6) 


and 
dW =rdz. 


Initial Orbits Circular 


When the orbit of the container is_ circular, 
de = ~/dh? + dl?, and dw = 0, from the definitions of 
h and 1. Consider a cylindrical shell, travelling on a 
circular orbit, whose (longitudinal) spin axis at the 
moment of explosion is yawed to starboard by a 


A 
Ww 


Fic. 2. Nomenclature and conventions 


degrees, and pitched upwards by 8. Here we have: 
VR = —pd¢[cos a cos ¢ + sin a sin ¢ sin B], 
v.s 
and 


V-w 


pdlsin a cos ¢ — cos a sin ¢ sin BI, 


I 


pdolcos 8 sin @ + } sin 2a sin 28 cos 4]. 

The eference direction ¢ = 0 lies in the plane formed 
by W and ® when the cloud is released. It is perpen- 
dicular to d, and points away from W, as in Fig. 2. 


Denoting by K the ratio of incremental speed to circular 
speed, 


K= "8 
na 


we suppose K* « K. Then, from Eq. (1), since uw = 0; 


ly =e 
5° 
dQ = 0, 
dt = K[cos 6 sin ¢@ + ¥ sin 2a sin 26 cos 4g], 
da = 2aKk|sin a cos ¢ — cos a sin ¢ sin BI, 
dh = K[cos a cos @ + sin a sin ¢ sin B], 
and 


dl = 2K sin a cos ¢ — cos a sin ¢ sin B]; 
and, consequently, 
2 


de’ = K*[cos a cos@ + sina sin ¢ sin 6] 


+ 4K"[sin a cos ¢ — cos a sin ¢ sin B). 


By Eq. (2), 

dx = K{[cos 6B sin @ + § sin 2a@ sin 28 cos ¢] sin u 
and 

ae = Uk 


From Eq. (3), 


dr = da — acos (u — w)de + 3(K’) ---, 


- =] dh 
@ = tan eT 


Substituting the proper expressions for da and de, 
we find that dr becomes: 


where 


dr = dall — sec w cos (u — w)]. 
Finally, we have from Eq. (6), to first order in K: 
dR = 2aK{sin a cos ¢ — cosa sin ¢ sin f] 
-[1 — sec w cos (nT — w)], 
dS = aK[sin a cos ¢ — cosa sin ¢ sin B} 


-[—3nT + 4 sec w sin (nT —w ) + 4 tan ol, (7) 
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and 
dW = ak{[cosB sing 


+ 3 sin 2a sin 26 cos ¢] sin nT’. 


Since a and 6 are independent and arbitrary, it follows 
that the secular term in dS vanishes when and only 
when a= 6 = 0. 


Some Examples 

By setting the eccentricity of the shell’s orbit equal 
to zero, we are free to examine several interesting ex- 
amples unencumbered by complex initial orbit pa- 
rameter effects. Three special cases illustrate the 
physical features of clouds which result from initial 
spin alignments purely in the S, W, and R directions. 

Case I. The spin axis is aligned with the velocity 
vector when the cloud is released. Then a = 6 = 0, 
w = +7/2, and Kq. (7) becomes: 


dR = —aK cos ¢ sin nT, 


ds 2aK cos ¢{1 — cos nT’, 


and 

dW = aK sn ¢ sin nT. 
Squaring and adding terms, we find 

(dR)? + (dW)? = (aK sin nT)’ 
and 


a = —tan ¢. 

Each particle harmonically traces a straight line in 
the R, W plane (Fig. 3), and the cloud itself becomes a 
radially pulsating cylinder, slowly oscillating about its 


A 


nT =37/g PS | 
nTt-=n/2 } 
— e | 


— 


ai 


i x 
\ =] aa i a a= . 
aI ae veal 

nl-r 


nt =79r/g 
nT=37/4 
nT=57/g 


Fia. 3. Pellet cloud history. Spin axis aimed forward at T = 0 — 


longitudinal axis, with its end-faces tilting sinusoidally 
with time (Fig. 3). 

Qualitatively speaking, the properties of this cloud 
are governed by the fact that the incremental velocity 
vector applied to each pellet has components in the 
W and R directions only. The contribution along W 
merely changes the inclination of the particle orbit, 
while the R-impulse modifies the eccentricity, but not 
the period, since an impulse applied normal to the 


T= 
ce NT=377/2 


nTl=27 


nT =47 


Pia. 4. Pellet cloud history. Spin axis aimed upward at 7 = 0 
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velocity vector leaves the energy, and therefore the 
period, unchanged. The combined perturbations are 
hus strictly periodic. 

Particles for which the R-impulse is directed inwards 
(ie., — 1/2 < @ < +7/2) find themselves at the 
latus rectum of the new orbit, travelling towards 
perigee; these move somewhat faster than nominally 
at first. Particles subjected to an outward-directed 
R-impulse (1/2 < ¢ < 32/2) move more slowly at 
first. Hence the exhibited periodic variation in the 
end-face slope of the cloud. 

4 Case IIT. When the spin axis is aimed northwards 
at T = 0, then a = 0, 8B = 7/2, and tan w = —} cot ¢. 
Again from Eq. (7): 


“dR = aK[2 sin ¢ (cosnT — 1) — cos¢ sin nT), 


dS = aK[— sing sin nT — 2 cos¢ (cos nT’ — 1) 

i + 3aK sin (nT — sin nT), 
and 

naW = (). 


Since dW = 0, there are no pellet plane changes, 
and the cloud cannot expand in height. This time, 
however, the spin motion also adds incremental speeds 
in the tangential direction. The energies are altered, 

and the cloud loses its property of periodicity. For 
particles located at ¢ = 7/2 and @¢ = —7/2, the 
motions are purely forward, and the pellets enter into 
their new orbits at apogee and perigee, respectively. 

Except for the secular term 3aK sin ¢ (nT — sin nT) 

appearing in dS, the figure traced by each pellet 
would be an ellipse with axes of length 2aKk, 4aK; the 


major axis of each ellipse is an extension of Ry, the 
particle’s initial radial position from the spin axis 
(Fig. 4). 

Ordinarily, therefore, the cloud shape would be not 
unlike the pulsating cylinder of Case I, but the secular 
term in this case superimposes a relative drift upon the 
motion of each pellet, so that the cylinder eventually 
dissipates itself into a smear along the orbit. Figure 4 
traces this development, starting from 7’ = 0. 

Case III. For a radially (outward) directed spin 
axis, a = 1/2, 8 = 0, and w = 0, or z. Here, 


dk = 2aK cos ¢[1 — cos nT], 

dS = aK cos¢sinnT — 3aK cos ¢[nT — sin nT}, 
and 

dW = aK sing sin nT. 


By analysis similar to the above, the cloud shape here 
exhibits the cylindrically pulsating, oscillatory be- 
havior discussed in Case I, coupled with a rectilinear 
drift of the type mentioned in Case II. 
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Estimation of Longitudinal Variations in 
the Earth’s Gravitational Field From 
Minitrack Observations 


W. M. Kaula? 


Abstract 


A term in the Earth’s gravitational field 
Vam = kM (Jnm COS MA + Kym Sin MX) (*) Pr» (sin ¢) 
Fs 


will cause a satellite to vary in its orbit with a frequency 
(0 or 1) + mQ — 6), where &, 2, and @ are the frequencies 
of rotation of the perigee, the node, and the earth respectively. 
Estimates of the magnitudes of the Jinn , Kin based on ter- 


1 Submitted March 27, 1961. 
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restrial gravimetry indicate that some of these orbital varia- 
tions will be of about +100 meters in amplitude. 

Investigation of other causes of residuals in satellite obser- 
vations—other orbital variations, electronic environmental 
and instrumental effects, orbital computation methods, geo- 
detic datum error—indicates that none of them will give rise 
to exactly the same spectrum of variations in the residuals as 
the gravitational harmonics. Hence the gravitational terms 
should be determinable from observations of adequate distri- 
bution and number. 

In the case of Minitrack observations of a single satellite, 
limitation of observations to high altitudes on the meridian 
makes it difficult to distinguish one gravitational term from 
another except through the difference in phase shift of the 
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spectrum for different stations. Hence the attempt to estimate 
the gravitational terms must be made simultaneously with 
the date from all stations. The results obtained from 3077 
observations of the Vanguard I satellite are as follows: 


Tag 300.69 212) 10 


Ja, = — 1.02 (+0.35) X 10~, 

Ke =i th ey Lore 
and 

Ku = + 0.30 (+0.35) x 107°. 


The principal defect in these results is probably systematic 
error in the orientations of the Minitrack antennas. 
Introduction 


The Earth’s gravitational potential may be repre- 
sented as a sum of spherical harmonics: 


V=- ae f = dls C) P.»(sin $) 


=I Cost hy Sin mo(4 ) Pan(sing) | 


n,m 


(1) 


in which kM = gravitational constant X earth’s mass, 
r = radius vector, ¢ = latitude, \ = longitude, R = 
Karth’s radius, Prm(sing) is a Legendre Associated 
Function, and the J’s and K’s are independent coeffi- 
cients. Jno is usually written as J, , and Kx» = O since 
sin 0 = 0. 

Jo = 1.0823 (+.0003) X 10°, the term due to the 
Earth’s oblateness, has been estimated now many times 
from secular changes of satellite orbits. Estimates have 
also been made of J,9 for n = 3 through 6, that of 
Jo = —2.3 (0.1) < 10° in particular being well 
confirmed (Refs. 1-4). 

However, in addition to these five J,,0’s which have 
been estimated, there exist 38 other coefficients Jnm , 
Kum, m # 0, S 6, which express longitudinal varia- 
tions in the gravitational field, and which are of just as 
great an interest to geophysicists. Existing terrestrial 
gravimetry indicates that the order of magnitude of 
these Jim , Kn», Will be the same as that already deter- 
mined for the J, : about 10 ° for normalized spherical 
harmonics. The terrestrial gravimetry is not distributed 
well enough, however, to give a reliable determination 
of any individual harmonic J,,, or Ky». Hence our 
purpose is to investigate whether it is feasible to esti- 
mate any of the longitudinal harmonics from satellite 
motions, and, if so, how. 

The expression in Eq. (1) of the gravitational poten- 
tial can be transformed from r, ¢, \ coordinates to 
osculating elements’ (a, e, 7, 2, w, MW) coordinates by 
some more or less tedious algebra (Ref. 5), ie., 


Z i A 2 
V= [2 > Goq(e) COs qM — ice (“) 
q=—0 i} 


a 


4 = Fop(2) ay Google) Sxopq(, M) 


p=0 q=—0 


2% al = ja 


a 


yd thay 


m,p=0 q=— 02 


Sanne (Cy Ils Gs 9) | : 
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where 
Ve (n—m) even 3 ) 
Samta =F ee cost (n aD | 
+ (n — 2p + q)M + m(Q — 6)} (3)) 


(n—m) even 
=|. ea sin{(n — 2p)w 
Dam (n—m) odd 
+ (n — 2p + g)M + m(@ — 8)} 


and a = semi major axis, e = eccentricity, 7 = inclina- 
tion, w = argument of perigee, J = mean anomaly,, 
and 2 = longitude of the node of the osculating orbit ;; 
6 is Greenwich sidereal time; and Frm,(7) and Gnpa(e)) 
are polynomials in the inclination and eccentricity, 
respectively. | 

The first term in V causes the unperturbed elliptical 
orbit. The most prominent effects of the /2o term are} 
the secular changes in the node, 2, and in perigee, @.. 
Since Jo is of the same order of magnitude as the Jnm i 
Kim, any attempt to deduce the latter from their first | 
order effects on an orbit must measure such effects with) 
reference to a computed orbit which takes into account: 
Joo at least to the second order, i.e., terms containing: 
J3) . In the present study, we used as such a reference: 
orbit the Hansen-type numerical general orbit de-. 
veloped by Musen (Ref. 6). 

The first order effects of the higher order coefiientl 
Jam, Knm on the orbital elements can then be deduced. 
by using the corresponding terms in V (with change of 
sign) as the disturbing function in the Lagrangian 
equations for variation of elements (Ref. 7) and inte- 
grating with respect to time. We then obtain, for exam- 
ple, for the perturbation of the node by a particular 
DALES im ies 


1 dF imp 


AQ aS 
Mi Ale 


esinipa au 


(4) 
Gave | Sate dt, 


where M is the mean motion. On performing the inte- 
gration with respect to time, | Srmpq dt, there appears 
in the denominator a factor 

[(n — 2p)a + (n — 2p + q)M + m(Q — 6)). 


If we take the unit of time such that 7 = 1 then, for 
example, for Vanguard I @ = 0.00181, 2 = — 0.00124, 
and 6 = 0.0932. Since these terms appear in the de- 
nominator, setting n — 2p + q = 0 amplifies the effect 
by 10/m; if further m = 0, the effect is either amplified 
by about 500 (m odd) or it is a secular perturbation 
(n even). As has been mentioned, this “cream” of small 
n,m = 0,n — 2p + q = 0 has already been skimmed 
off, so we must turn our attention to the pee of 
msmall,n — 2p + q = 0: to begin with m = 
Crna} is of O(e''), and for n odd, q # Oforn — 2p + 
= 0: hence the longitudinal terms of greatest effect 
abate have n even. Jo, , Ko; cannot exist, since the 


Earth must rotate about an axis of maximum moment 
Wf inertia; hence the coefficients most likely to be effec- 
tive are J , Ka . Taking estimates of the magnitude of 
these coefficients based on autocovariance analysis of 
terrestrial gravimetry (Ref. 8), we obtain as estimates 
of the magnitude of effect of Ja, Ka on Vanguard I 
about +200 meters. Most of this effect is in the perigee 
argument and mean anomaly—along the orbit; but in 
this respect we benefit from sole advantage of the short 
period over the long period analysis: the drag effects of 


Other effects estimated from the terrestrial statistical 

data are Jo, Ko: +100 m.; Jn, Ku: +100 m; 
PYe2, Ke : +100 m; Je, Ke : +50 m. (Ref. 5). 
_ From expressions such as Eq. (4) we can set up a 
6 X k matrix for the partial derivatives of the osculating 
elements with respect to k of the gravitational co- 
Meiticients Jnm, Knm : 


d( Elose) 
Oda 5) e=) ; 


a By applying the appropriate rotations, we can obtain 

the position vector X of the satellite in a coordinate 

| system referred to the Earth’s axis and the vernal 

equinox from the vector Q referred to the orbit: 

ex = F;(—2)R,(—7)R,(—2)Q, 

b a(cos H — e) ) (6) 

Q= SVE RB 
0 


Cre = (5) 


The index indicates the axis for a rotation matrix 
and the argument indicates the angle rotated, with + 
“for a counterclockwise rotation. By differentating 

Eq. (6) we obtain the partial derivatives of the com- 
ponents of X with respect to the osculating elements: 


-_ 


Bile OCs, Ye) a 
Cxz = hla : (7) 
Rotating angle @ about the z axis from the vernal 
equinox to the Greenwich meridian and subtracting 
the rectangular geodetic coordinates U» of an observing 
station gives the satellite position in a geodetic system 
referred to the station. Applying further rotations 
gives the satellite position referred to axes peculiar 
to an instrument at the station; for example for a 
radio interferometric system with east-west and 
north-south baselines: 


ar, = R, (Z) x. (G o- 2) R;3(d) [R3(T).X — Uo 


= Rig|Rs(T)X =a Ud]. 


To obtain finally the actual observed quantities one 
more matrix N must be applied to L; for the direction 
cosines 1, m with respect to the E-W, N-S baselines, 


(8) 


Nix ae _ (9) 


where p is the range from the station to the satellite. 
Differentiating Eq. (9) with respect to the elements 
of L and combining it with Eqs. (5), (7), and (8) 
we obtain an observation equation relating corrections 
to the direction cosines 1, m to corrections to the station 
position and to the gravitational coefficients: 


wiuns(E) fg #) 


| ReloCxal0, z, oM)Caay eb 


-aw]-a{)) 
Mm) oBs 
thee 
M) oBs M) come’ 


where we have indicated the angular arguments in Cxz . 

The foregoing analysis was originally made in 
anticipation of a proposed geodetic satellite to be 
observed several times a day with effective accuracies 
on the order of +30 meters, from which several gravi- 
tational coefficients should be deducible. What we 
actually have in hand are Minitrack observations of 
about -+-1 mil residuals—or an effective accuracy on 
the order of +1000 meters. About 6000 of these ob- 
servations on Vanguard I have been accumulated over 
almost three years, so if the +1000 m. errors were 
perfectly random, we could get good estimates of 
several of the Jnm, Knm - 

However, the errors are not, of course, random, so 
we must inquire into what is the probable spectrum of 
the errors and how much does this spectrum overlap 
or alias the spectrum of effects of the gravitational 
coefficients Jim, AKnm we are attempting to estimate. 
We can sort the possible sources of error into six 
classes according to the sequence of events from the 
original motions of the satellite in its orbit to our 
estimation of the Jam, Knm: (1) Variations of the 
satellite in the orbit due to physical effects not fully 
allowed for, such as drag, radiation pressure, tidal 
effects, and gravitational coefficients not included; 
(2) Variations of the satellite-to-station path of the 
radio signal from that assumed, due to irregularities of 
ionospheric and tropospheric refraction; (3) Geo- 
metrical limitations on the time and place of observa- 
tion; (4) Tracking station position or orientation 
error; (5) Variations in the internal workings of the 
Minitrack instrumentation; and (6) Variations of the 
computed direction cosines due to the procedures 
used in the orbit computation, such as switching to a 
new set of reference elements each week, neglecting 
rotation of the orbit in differential corrections, etc. 
Each of these possible sources is most easily visualized 
in a particular coordinate system, to which our obser- 
vation equation (10) or vector equation (8) can be 
transformed by application of the appropriate rotation 
matrices. 

Classes (1) and (6) are most conveniently examined 
in terms of the orbital elements. As we have discussed, 


(10) 
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the orbital variations on which estimation of Jum , 
Um depends are primarily those of frequencies m(Q — 6) 
for n even; +20 + m(Q — 6) frequencies have co- 
efficients of O(e’). In the particular case of Vanguard 
I and nm = 41, the primary frequency is 1.014 cycles/ 
day, while the side frequencies are 0.976 and 1.052 
cycles/day. Exactly the same frequencies will arise 
from nm = 61, 81, ete. Although the amplitudes of 
effects on different orbital elements of one satellite 
will differ, a simultaneous least squares solution for all 
these coefficients from one satellite is impractical 
because the non-uniform distribution of observations 
makes the separation of these effects weak. Slightly 
different frequencies of 0.995 and 1.033 cycles/day 
will arise from nm = 31, 51, ete., with coefficients of 
O(e). 1/(1.014-0.995) = 53 days, so the effects should 
be easily separable with the available data. 

There are a host of gravitationally caused variations 
in the orbit of frequency M = 10.7 cycles/day (in- 
creasing 0.015 cycles/day/year), 21.4 cycles/day, 
32.1 cycles/day, ete. with total amplitude less than 
+100 meters (Ref. 5). The frequency of observations 
by the Minitrack averages about 6/day, giving a 
“folding” frequency of 3/day and aliasing frequencies 
for 1 cycle/day of 5, 7, 11, 13, 17, 19, ete. cycles/day 
(Ref. 9). None of these frequencies are close enough 
to the 10.7 cycles/day, etc. to be a likely cause of 
difficulty for the long record of data available, how- 
ever, the non-uniform distribution of observations 
may enhance the distortion of estimates of coefficients 
of lower frequency variations. 

Luni-solar attraction, radiation pressure, drag, and 
Earth tides do not appear to give rise to any marked 
daily or semi-daily periods in the orbit, (unless there 
is a pronounced longitudinal variation in atmospheric 
density due to the funneling effect of the geomagnetic 
field on solar charged particles). The dominant effect 
of drag is best characterized as a “‘noise’’ of continuous 
spectrum, from which the discrete spectrum caused 
by the gravitational terms should be resolved by a 
sufficiently long record. 

In class (6) error might possibly arise because the 
reference orbits are computed for seven-day periods. 
The orbit determination finds a secular drag coefficient 
and parameters of a rotating Hansen ellipse such 
that, in conjunction with specified values of kM, Ju, 
and J4 variations of the satellite position are most 
closely accounted for. The most prominent of these 
variations are expressed as secular changes in the 
cyche variables M, w, Q; if the net changes in these 
variables over the one-week period were in part due 
to a Jum, Knm which we are estimating, then the orbit 
determination procedure will contaminate estimates 
made from residuals with respect to the orbit. This 
contamination would be a maximum for Jam, Knm 
causing variations of frequency (22 + 1)/2W, where 
I is an integer and W is the duration covered by the 
orbital determination, and would be zero for Jnm 
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causing variations of frequency J/W. (Ref. 10). Fo 


precisely seven-day periods without overlap contami} 
| 


nation peaks should exist for 2.071, 1.927, 1.071, 0.924 
eycles/day: all of which are fairly distant from the 
frequencies of primary interest: 2.028 and 1.014 


cycles/day. | 

Errors in classes (2), (4) and (5) the electronic 
environmental effects position and orientation error 
and instrumental effects, are most easily visualized as 
variations in the observed direction cosines 1, m. 
should expect that the most prominent peak in those 
error spectra will be close to 1.000 cycles/day, because 
of variations in the ionospheric and _ tropospher ic 
refraction, the geometrical limitations on time of ob- 


servation, class (3), and the temperature and humidity 


Minitrack stations. Since most of these variations a er 
all rather irregular, we should expect this peak to bes 
rather broad. The frequency in which we are mos i 
interested is 1.014 cycles/day in terms of the orbital] 
elements; to determine whether it might be contami- 
nated by the refraction effects, however, we must} 
first apply the various rotation matrices in Eq. (10)} 
tod = Cea pm 5 tinewte 7 

Since Eq. (10) is rather complicated, we shall make 
this transformation only to the first approxima-- 
tion—i.e., drop all terms with coefficients including: 
the eccentricity or powers thereof, and make a further’ 
abbreviation in the form of a notation: 


(q1 9243 2 -) = Tso {q@it+e} 


sin 


(11), 


where K = [Jiik; and; = []in (4a +@ + =f 
+q,). Then 


dE” = {(mQ — mé@), (mQ — me), (m2 — mé), 0, 


(mQ — mé), (mQ — mé)} 
and 
Q” = {(), (M), 0}, ete. 


In the various rotations applied in Eq. (10), 
Ri(—2), Rs), R3(r/2 — @) and R3(/2) have constant 
arguments, and therefore constitute phase shifts which 
do not introduce any new frequencies. As it turns out, 
these phase shifts are important to making a solution 
possible, but for the present we shall ignore them 


sO 


that transforming from dE to {i is the same as to 


i 
re so far as frequencies are concerned. Carrying 


out the various multiplications, we get, substituting 
S = mQ — mé: 


l 1 
ae comm = pte = 1 (0) + (Ga) + (wM) 
(12) 


+ (68002AM M) + (wo) + a 


where 
th = dl, = (@SMQ) + (6SMw) + (6SMaQ) 
+ (@8SwQ2MM) + (@Sw02) + (SM) + (Sel) 
+ (Sw) + (SoMM) + O(e). 


Before carrying out the tedious multiplications, we 
can simplify by applying the various conditions of 
observation, which we have called class (3): 

(1) For a given Minitrack station, the satellite is 
only observed crossing the meridian plane; i.e. for the 
right ascension, a, of the satellite, 


Ca ake (13) 


cot (w + f) = cot (A+ @—)cos7z. (14) 


For Vanguard I, neglecting terms of O(e) = 0.19 for 
Vanguard I means setting M = f and cosi = cos 


34° = 1, whence 
a M=2ir+60-2—w0+ Ole) (lis) 


: 
(2) At a given Minitrack station, the satellite is 
observed only when it is within about 22° of the zenith; 
so that 


| $ (Satellite) 
— ¢ (Station)| < 2.1° (16) 


at perigee, 


and at apogee, | ¢ (Satellite) 
— ¢ (Station)| < 8.8° 


¢ (Satellite) = ¢ (Station) + O(e) 
and 


f=sint(S2¢) — w+ O(c) (17) 
sin 7 
so 

M = —# + O(e). (18) 
Combining Eqs. (18) and (15) gives: 


6= 7+ O(e) (19) 


and 
S=0+4 0(e). (20) 


Equation (20) indicates that at a single station all 
eravitational perturbations due to terms containing 
argument md look alike, to the first order; if they are 
io be distinguished, it must be by terms of O(e) and 
higher. It indicates further that the solution for several 
of the Jnm, Knm Will depend on the differences in 
phase shift introduced for different stations by the 


rotations re( — ‘), R;(\), and that the solution for 


the gravitational coefficients should be made simul- 
taneously from several stations by least squares, 
rather than by harmonic analysis. Substituting Eqs. 


(18), (19), and (20) in (12), we get: 
: = (0) + deny ron) Cron ae 
+ (www) + O(e) 


= 2 20 + jw) ai O(e) (21) 
dL = (QQw) + (Qww) + (QQew) + (NQwww ) 
+ (QQw) + (w) + (ww) + 


cee 38 


+ (www) + O(e) = Xd (g2+hw) + O(e). 


The ; periods in Eq. (21) will appear for all effects 


which affect the direction cosines through the satellite 
range, which would include the datum error, dU, , 
and, to some extent, the ionospheric refraction error. 
The dL periods will be peculiar to even-degree longi- 
tudinal harmonics in the gravitational field; it is the 
periods that make these harmonics separable from 
other effects—in particular, from the electronic en- 
vironmental and instrumental effects, for which we 
can expect daily, monthly (change of refraction for- 
mula), possibly semi-annual (recalibration), and 
annual periods. 

In view of the foregoing considerations, a least 
squares solution was attempted using observation 
equations of the type of 10 simultaneously for the ten 
even-degree harmonics likely to have greatest effect: 
Jo, Ko, Jn, Kun, Je, Keo, Ja, Ka, Jeo, Keo, 
which probably include all even-degree terms having 
effects of +20 meters or more amplitude in daily and 
shorter periods. (J 3, , Us, might have a greater effect, 
as previously mentioned, but is is separable due to the 
extra w term in its argument.) However, the result 
indicated weak conditioning: in general, the combined 
effect of a set of similar coefficients (e.g., Jo, Ja, 
Je.) Was appreciably smaller than their separate 
effects, and the coefficients were unreasonably large 
compared to the statistical estimates from terrestrial 
gravimetry (Ref. 8). So solutions were made for J2 , 
Ky, Ja, Ku only, and results were obtained of rea- 
sonable magnitude: 

From obs. 1—-1128: 


Jo, = +0.27 (40.59) X 10°, 
Ko = +2.59 (40.61) X 10° 
Jn = —1.67 (40.20) x 10°, 
and 
Ka = —0.09 (+0.23) xX 10°. 
From Obs. 1129-2241: 


Jomo Ao (e051 ele 
Ko = +0.85 (40.58) X 10° 
ee ONE 22) Xa Or, 
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and 
Ka = +.53 (40.23) X 10°. 
From Obs. 2242-3077: 
Jog == 0.99.60, 75 eX lan 
Kor= 08! (067 ane 


Jig a 0865 0. 26)E xX 108. 
and 
Kua = +0.52 (40.26) X 10°. 


The uncertainties stated are those obtained in the 
conventional manner from the least squares solutions. 
The changes in results from one solution to another 
are disappointingly large. These observations covered 
525 days. The Minitrack stations are recalibrated 
about once each six months; changes in antenna 
orientation on recalibration have been as great as 
4 X 10“, so the foregoing changes in results could be 
due to variations in the effective orientations of the 
Minitrack antennas. The uncertainties for a combined 
solution are therefore better taken as based on the 
discrepancies between separate sets: 


(SEO ali 
Ke = +117 C212) « 10% 
Ja = —1,02-(+025).x 10 


and 
Ke = -+0.30 (£0.35) < 107°. 


Provided the antenna orientations remain stable, 
Minitrack data should be adequate to determine the 
leading coefficients of the gravitational field. In any 
case, whatever, the type of observations used, more 
analysis in the direction started in this paper is needed 
to extract all the information existing in satellite 
observations. In particular, more investigation is 


needed of the interaction between the gravitationally 
caused variations in satellite orbits and the position o¥ 
tracking stations, and of the effect thereon of the non 
uniform distribution of observations due to tk 
gesonetrical limitations of the observing system an 
orbit used. 
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The rigid body assumption in the solution of motion of 
bodies m space, though good as a first approximation, has 
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energy dissipation, one has to extend the analysis and con- 
sider elastic deformations. Using expressions for gyroscopic 
forces derived from the rigid body assumption, an elastic solu- 
tion for the vibration of a disk is obtained. Analytical expres- 
sions for the deflection, bending moment in the radial diree- 
tion, and vertical force were derived and plotted as a function 
of the radial distance. 


I. Introduction 


The motion in space of a moment-free rigid body of 
revolution, with principal mass moments of inertia 


A, A, C, consists of steady precession of the symmetry 
ixis about the angular momentum vector, h, fixed 
in space. The angle between the symmetry, or spin 
uxis, and the direction of the vector h, is called at- 
itude angle and is a constant. A real body, however, 
is elastic. A vibrating elastic body undergoes stresses 
and deformations resulting in internal energy dissipa- 
tion, which, in turn, causes a change in the attitude 
angle, 6. The assumption is made here, that the 
elastic deformations are small relative to the unde- 
formed geometry of the body. Furthermore, it is 
assumed that the rate of change, 6, or nutation velocity, 
is small relative to the precessional velocity, ~, and 
“spin velocity, ¢. The forces causing the elastic vibra- 
tions are the gyroscopic or inertial forces. 
_ The problem to be discussed here treats the case of 
‘two circular elastic disks connected by a rigid shaft 
(Fig. 1). 
3 The solution of a circular plate having the inner 
circular edge clamped, the outer one free, and sub- 
jected to gyroscopic forces is to be found. In order to 
‘determine the gyroscopic forces, an expression for the 
acceleration at any point, P, must be derived. 

The acceleration is expressed in terms of the initial 
angular velocity, wo , and the attitude angle, 6, treated 
as a parameter. 


Il. The Acceleration Expression 

Since only the rotational motion is of interest and 
there is no relative motion in rigid bodies, the ac- 
celeration of any point, P, is given by: 

ap =o XoXR+0XR. (1) 

where w and are the angular velocity and acceleration 
of the body respectively, and R is the radial distance 
fo P. 


In terms of components along the body axes, 1, 2, 
3, the radius vector R is given by: 


R = (rcosa)i + (r sin a)j + 2k. (2) 


>| 


p - precession 


Fra. 1. The body and its angular velocities 


Fic. 2. Coordinates of point P 


In Fig. 2, the system XYZ is fixed in space and 
the angular momentum vector, h, is taken along the 
Z axis. The xyZ system is obtained by rotation in the 
XY plane by an angle, y. System zy’z results from 
rotation of system «yZ about the node axis, by an 
angle 6. System 1, 2, 3 is attached to the body and 
rotates with respect to system xy’z with angular velocity 
¢ while the direction of axis 3 coincides with the one of 
z at all times. It follows from Figs. 1 and 2 that the 
angular velocity components along the body axes, are: 


w = y sin 6 sin ¢, 


w. = W sin 6 cos ¢, a 


and 
wo = d+ wv cos 0. 
EKuler’s equations for rotational motion in space 


can be written in vector form, 


d / 
M = dt (h), (4) 


in which M is the external moment vector and h the 
angular momentum vector given by 


h = A (ai oa wo ) + Cousk. (5) 


In the above, A and C are the mass moments of inertia 
about the pitch and spin axes, respectively. 

It can be shown from Hqs. (4) that, if there are no 
external moments applied: 


w3 = const, (6) 
6 = const, Ge. 
Co 
mewn ee 8 
v (A — C) cos @’ (8) 
and 
; (C W3 
=e 9) 
u A cos 0 (9. 


In reality, 6 is not zero but very small compared to 
WY or ¢. With this in mind, the angular acceleration 
components become: 


a, & vd sin 0 cos ¢, 
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»» & —y¢ sin 4 sin ¢, 

@ & —yo sin 6 sin ¢ (10) 
and 

a; = 0. 


The precession and spin velocities ¥ and ¢ can be 
expressed in terms of the initial angular velocity, w . 


It was shown (Reference 1) that if _ > 1 then 6 


> 0. The case where the spin axis initially coincided 
with the direction of the angular momentum vector, h, 
is considered here. This implies that, initially, @ ~ 0 
and the magnitude of the constant angular momentum 
vector is | h| = Ca. 

From Fig. 3, the following relations can be written: 


li, Shy oosa ee Cun cos 6 = Cus , 
and (11) 
Inn = | h| sin 6 = Cup sin 6 = A/a? + a. 
One concludes that: 
3 = wo COs 0 (12) 

and since 

Vor + o? = sin 8, 
it follows that: 


eG 
y Pe a wo : (13) 
Substituting Eq. (13) into (8), one obtains 
C 
o= (1 — a Cos 6. (14) 


Expressions (13) and (14) introduced into (3) and 
(10) yield the angular velocity and _ acceleration 
components in terms of wo. Finally, this leads to the 
acceleration vector, ap , if Eq. (1) is used. 


2 
aes wi{] - (4) sin’ écos¢cos (a+ ¢) 


2 
—rcos 6cosa +z @) sin 6 cos @sin 6| i 
ONG ra 
+ Ir 4) sin Osin dcos(a+ ¢) 
ne (15) 
—rcos sina +z ©) sin 6 cos 6 cos 6 |i 


abe E o (2 me & sin Ocos @sin (a + ¢) 
A A 


—2Zz (G). sin’ | k} 3 


IIT. Elastic Solution of the Plate 

The exciting forces on the system are the inertia or 
D’Alembert forces. The deformed shape is assumed 
very close to the undeformed one, so that the real 
accelerations are approximated by the rigid body 
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Fig. 3. Angular momentum vector components % 
accelerations (given by Eq. (15)). Since only the} 
component causing bending is of interest, only the} 
component in the direction of the unit vector k ish 
retained. From this, only the portion causing vibratory 
motion is of further interest. Note that the oscillatory 
terms are the ones containing the angle ¢ = ¢t, where 
é is the spin velocity and ¢ denotes the time. Hence? 


ap Si (2 - ) sin @cos@sin (a+ ¢t)k. (16) 


Denoting by qg the distributed force per unit volume 
and by p the mass density one writes: 


q = pap = ane (2 — o) sin 6cos @sin(a + ¢t) 


(17) 
= frsin (a+ ¢t) 
where: 


, oC ae 
fF = poo a (2 _ “) sin 6 cos 6. (18) 


i 


As far as the solution of the plate is concerned, f 
represents a constant. 

The distribution of the gyroscopic forces is shown 
in Fig. 4. 

Taking into account the direction of the exciting 
forces, the plate equation is written: 


Vi w+ a-rsin (a+ 6¢) +5 = (19) 
in which Dz, is a modified flexural rigidity given by 
D Eh 
D 7 = .- = —___,, 92 
"ate Om 


Polar coordinates should be used, due to the shape 
of the boundaries. The variables in Eq. (19) can be 
separated by the following transformation: 


w = W(r) sin (a + ¢t), (21) 
where W(r) is a function of r alone. It will prove 


convenient to introduce a new parameter, y, given by: 


og = 7 (22) 


Fia. 4. Gyroscopic force distribution 


Introducing Eqs. (21) and (22) into (19), one ob- 


-; 
_ tains an equation for W(r) alone, which when it is 
solved, gives: 


Wr) = X (1,1) Ai(yr) + X(2,1) Vilar) 
: 


© +X, hOr) + XG1)K(r) + iD 


(23) 


r. 


In the above, J:(yr) and Y,(yr) are Bessel functions 
of the first order and first and second kind, respectively. 
—I,(yr) and K,(yr) are modified, or hyperbolic Bessel 
functions of the first order and first and second kind 
“respectively. The coefficients, X(j, 1), have to be 
determined from the boundary conditions. Combining 
expressions (21) and (23), the solution of Eq. (19) 
is written: 


w=[ XU, DAG) +X QD NG") 
+ X(3,1)h(yr) + X(4,1) Kilyr) (24) 
Hi | Ves 
Sear t). 
msi.” sin (a+ ¢t) 
Solution (24) is subject to boundary conditions. 
Fig. 4 indicates that at the inner boundary, r = 6, 
where the plate is clamped, the deflection and the 


slope to the deflection curve in the radial direction 
are zero. 


w|a» = 0 (25) 
| 
Oe carers (26) 
or i) 
At the outer boundary, r = a, the plate is free, 


so that the bending moment in the radial direction 
and the vertical force on areas normal to r are zero 


aw low , 1 a) 
A, = =0 (27 
ee E ae (: or . r dee ind 2m) 


and 


(28) 


1 dw a4) ) = 0 
r Orda ~ r 0a/ |\r=—a 


The four boundary conditions, Eqs. (25), (26), 
(27) and (28), yield four simultaneous equations 
for the coefficients X(j, 1). The four equations are 
of the type: 


A(t, j)X(j, 1) = BC, 1), (29) 


where summation over j is implied. 
The boundary condition Eq. (25) leads to: 


ACT 1) = Jalxb), 
A(1, 2) = Yi(yb), 
A(1, 3) = h(yb), 


A(1, 4) = Ki(yb), (30) 
and 
ee 
B(1,1) = ape b. 


From boundary condition Eq. (26), one obtains: 
A(2, 1)'= ybJo(yb) — Jilyo), 
A(2,2) = ybYo(yb) — Yilyb), 
A(2,3) =—syol(yo) —itge)) 
A(2,4) = —ybKo(yb) — Ki(yb), G1) 


I 


and 


j 
B(2 =— : 
(2,1) = — 


| 


After more laborious manipulations, boundary 
condition Eq. (27) yields: 


AO el ve el ze) Jo(ya) 
ya 
2(1 — pv) 
st i a | Ji(ya), 
A(3,2) = eee Yo(ya) 
ya 
2(1 — v) 
aie pee Fe 1| Yi (ya), 
A(3,3) = ~via) Iy(ya) 
2(1 ) Sa 
= Eee ot 1] I(ya), 
A(3,4) = eres) Ko(ya) 
a 
i ES ee) ie | Kiya), 
(ya) 
and 
BG, 1) = 0. 


Finally, after simplifications, boundary condition 
Eq. (28) results in: 


A(4,1) = -1 a | J (ya) 
(ya)? 
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1 21 — a 
ar ya E -- ~ (ya)? Ji(ya), 
A(4,2) = -1 + a aE | Yo(ya) 
1 2(1 — zt 
4 na E == eat Yi(ya), 
Cl 7) 
A(4,3) =| 1 — ~~ | hlve) 
| (ya) | | . (33) 
1 2(1 — » 
~ ya [: Ga? | noe, 
A(4,4) = E wee Ko(ya) 
1 2 v) 
Fel? ap |B, 
and 
BIA) =O: 
The set of Eqs. (29) can be written in matrix form: 
A eS = By, (34) 


in which [A] is a square matrix, and {X} and {B} 
are column matrices. 

The solution of Eq. (384) is readily obtained if one 
premultiplies both sides by the inverse of matrix [A], 
Le., 


{X} = [A] {B}. (35) 


The above expression gives the coefficients, X(j, 1), 
which, when introduced into Eq. (24), yields the 
complete solution for the deflection, w. 

Making use of Eqs. (27) and (28), one can write 
equations for the bending moment, M, , and vertical 
force, V, : 


M, = yD Ea {-" 


2(1 — v) : 
+ Ga 1] 409} 


+ X(3,1) ee? faces (36) 
yr 


Be ae ae | hian)} 
Gray 


+ pee + 1| Kon) | 
(yr)? 
- sin(a + ¢t) 
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Attitude Angle 6 = 45° 


W (r) - in 


M, (r) - Ib in/in 


V(r) vs 6 

c 
= 0 
2 
50 + | 
> 100+ | 
> 2 

200 

0 4 6 8-10 12° 945 62) Ses 
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IV. Numerical Example 


The amplitudes of the functions w, M,, and V, 
were evaluated and plotted for the following values: 


Attitude angle 6 = 45° 

External radius a7=" 20.030 

Internal radius bv= 4 05n 

Thickness h =O) an 

Mass density p = 0.00073237 Ib in® 
sec” 


Poisson’s ratio vy =-0.30 
Young’s modulus E = 30.710" psi 
_ Initial angular ve- wy = 20 rad sec * 
locity 
Moments of inertia Chea = 0.1 
ratio 


The results are shown in Fig. 5. 


V. Conclusions 


In many dynamical problems, the rigid body solu- 
‘tion can be regarded as a first approximation only. In 
certain cases, such as with space vehicles, the rigid 
_ body analysis is not sufficient and one has to consider 
the elastic vibrations. In many cases, the acceleration 
given by Eq. (15) can be used for determining the 
inertial forces. The assumptions made in the present 
i analysis were the conventional small deflection assump- 
tions and small nutational velocity, 6. Both assumptions 
are justified except when resonance becomes important. 
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VII. Symbols 


Symbol Definition 
AG, = moment of inertia about pitch and spin axes re- 
spectively 
h = angular momentum vector 
6 = attitude angle measured from vector h to spin axis 
6 = nutation velocity 
> = spin velocity 
~y = precession velocity 
18 = arbitrary point 
ap = acceleration vector 
ra) = angular velocity vector 
r6) = angular acceleration vector 
R = radius vector 


1,2,3 = body axes 


i,j,k = unit vectors along body axes 

r,a,z = polar coordinates 

X YZ = fixed axes in space 

xyZ = coordinates system 

zy'z = coordinates system 

W1,9,3 = angular velocity components about body axes 

M = moment vector 

O1,2,3 = angular acceleration components about body axes 

hav. = angular momentum components 

| h | = angular momentum magnitude 

wo = initial angular velocity 

p = mass density 

q = distributed force 

w = plate deflection 

Wir) = amplitude of plate deflection 

D = plate flexural rigidity 

E = Young’s modulus 

h = plate thickness 

v = Poisson’s ratio 

a, b = outer and inner radius of plate 

M, = bending moment in radial direction 

V; = vertical force on plane with normal in the radial 
direction 
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